Enhanced binding and cold compression of nuclei due 
to admixture of antibaryons 
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Abstract 

We discuss the possibility of producing a new kind of nuclear system by putting a 
few antibaryons inside ordinary nuclei. The structure of such systems is calculated 
within the relativistic mean-field model assuming that the nucleon and antinucleon 
■ potentials are related by the G-parity transformation. The presence of antinucleons 

CS . leads to decreasing vector potential and increasing scalar potential for the nucleons. 

As a result, a strongly bound system of high density is formed. Due to the significant 
reduction of the available phase space the annihilation probability might be strongly 
suppressed in such systems. 



PACS: 25.45HL 27.20.+n, 21.10.Gv, 21.30.Fe 
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Presently it is widely accepted that the relativistic mean-field (RMF) model |IJ gives a 
good description of nuclear matter and finite nuclei || . Within this approach the nucleons 
are supposed to obey the Dirac equation coupled to mean meson fields. Large scalar and 
vector potentials, of the order of 300 MeV, are necessary to explain the strong spin-orbit 
splitting in nuclei. The most debated aspect of this model is related to the negative- 
energy states of the Dirac equation. In most applications these states are simply ignored 
(no-sea approximation) or "taken into account" via the non-linear and derivative terms of 
the scalar potential. On the other hand, explicit consideration of the Dirac sea combined 
with the G-parity arguments leads to such interesting conjectures as the existence of 
deeply-bound antinucleon states in nuclei || or even spontaneous production of nucleon- 
antinucleon pairs 0, f| . These predictions are based on the assumption that the relativistic 
description of the nucleon within the RMF model is a valid concept (for a discussion see 
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Refs. |6|, 0, ||). Unfortunately, the experimental information on the antinucleon effective 
potential in nuclei is obscured by the strong absorption caused by annihilation. As follows 
from the analysis of Ref . || , the real part of the antiproton effective potential might be 
as large as 200-300 MeV, with the uncertainty reaching 100% in the deep interior of the 
nucleus. 

Keeping in mind all possible limitations of the RMF approach, below we consider 
yet another interesting application of this model. Namely, we study properties of light 
nuclear systems containing a few real antibaryons. At first sight this may appear ridiculous 
because of the fast annihilation of antibaryons in the dense baryonic environment. But 
as our estimates show, due to a significant reduction of the available phase space for 
annihilation, the life time of such states might be long enough for their observation. In 
a certain sense, these states are analogous to the famous baryonium states in the NN 



system ||10|| , although their existence has never been unambiguously confirmed. 

It should be emphasized that previous discussions of the negative-energy states in 
nuclei were somewhat academic because their contributions to the source terms for the 
meson fields were ignored. In contrast, our primary goal here is to study properties 
of physical systems containing both real baryons and real antibaryons. There could be 
different theoretical schemes for solving this problem but we find the RMF model most 
suitable for this study. To our knowledge, up till now a self-consistent calculation of 
antinucleon states in nuclei has not been performed. Our calculations can be regarded as 
the first attempt to fill this gap. We consider two nuclear systems, namely 16 and 8 Be, 
and study the changes in their structure due to the presence of an antiproton. 



2 Theoretical framework 

Below we use the RMF model which previously has been successfully applied for describing 
ground-states of nuclei at and away from the /^-stability line. For nucleons, the scalar and 
vector potentials contribute with opposite signs in the central potential, while their sum 
enters in the spin-orbit potential. Due to G-parity, for antiprotons the vector potential 
changes sign and therefore both the scalar and the vector mesons generate attractive 
potentials. 

To estimate uncertainties of this approach we use three different parametrizations of 
the model, namely NL3 [|ll|], NL-Z2 [[12] and TM1 |13 |. Their parameters are found by 



fitting binding energies and observables related to formfactors of spherical nuclei from 16 O 
(not included in the TM1 fit) to Lead isotopes. In NL3, properties of symmetric nuclear 
matter have been included in the fit as well. The TM1 model, implementing a self- 
interaction of the cu-meson, gives a softer rise of the vector potential with density. This 
leads to smaller meson fields as compared to the NL3 and NL-Z2 parametrizations. In this 
paper we assume that the antiproton interactions are fully determined by the G-parity 
transformation. Following Mao et al. [O], we solve the effective Schrodinger equations for 
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Figure 1: The left panel represents the sum of proton and neutron densities as function 
of nuclear radius for 16 without (top) and with an antiproton (denoted by AP). The 
left and right parts of the upper middle panel show separately the proton and neutron 
densities, the lower part of this panel displays the antiproton density. The right panel 
shows the scalar (negative) and vector (positive) parts of the nucleon potential. Small 
contributions shown in the lower row correspond to the isovector (p-meson) part. 



both the nucleons and the antiprotons. Although we neglect the Dirac sea polarization, 
we take into account explicitly the contribution of the antibaryon into the scalar and 
vector densities. For protons and neutrons we include pairing correlations within the BCS 
model with a 5-force (volume pairing) |15|. Calculations are done within the blocking 
approximation [|l(J for the antiproton, and assuming the time-reversal invariance of the 
nuclear ground-state. The energy of the system is found by using the damped gradient 
iteration method . The coupled set of equations for nucleons, antinucleons and mesons 
is solved iteratively and self-consistently. The numerical code employs axial and reflection 
symmetry, allowing for axially symmetric deformations of the system. 



3 Structure of nuclei containing antiprotons 

As an example, we consider the nucleus 16 with one antiproton in the lowest bound 
state. This nucleus is the lightest nucleus for which the mean-field approximation is 
acceptable, and it is included into the fit of the effective forces NL3 and NL-Z2. The 
antiproton state is assumed to be in the sl/2 + state. The antiproton contributes with the 
same sign as nucleons to the scalar density, but with opposite sign to the vector density. 
This leads to an overall increase of attraction and decrease of repulsion for all nucleons. 
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Figure 2: Proton (left), neutron (middle) and antiproton (right) energy levels for the 
nucleus 16 O with one antiproton and without it (rightmost columns in left and middle 
panel). 



The antiproton becomes very deeply bound in the sl/2 + state. To maximize attraction, 
protons and neutrons move to the center of the nucleus, where the antiproton has its 
largest occupation probability. This leads to a cold compression of the nucleus to a high 
density. 

Figure [1| shows the densities and potentials for 16 with and without the antiproton. 
For normal 16 all RMF parametrizations considered produce very similar results. The 
presence of an antiproton dramatically changes the structure of the nucleus. The sum of 
proton and neutron densities reaches a maximum value of (2 — 4) po > where po — 0.15 fm -3 
is the normal nuclear density, depending on the parametrization. The largest compression 
is predicted by the TM1 model. This follows from the fact that this parametrization gives 
the softest equation of state as compared to other forces considered here. 

According to our calculations, the difference between proton and neutron densities is 
quite large, which leads to an increase in symmetry energy. The reason is that protons, 
though they feel additional Coulomb attraction to the antiproton, repel each other. As 
a consequence, neutrons are concentrated closer to the center than protons and the sym- 
metry energy increases. As compared to the case without the antiproton, absolute values 
of vector and scalar potentials increase in the central region of the nucleus. This leads 
to an enormous drop of the effective baryon mass near the nuclear center, which should 
strongly suppress local rates of annihilation (see below). The effective mass for TM1 even 
becomes negative at r < 1 fm. 

Since nucleons feel a deeper potential as compared to the nucleus without the antipro- 
ton, their binding energy increases too. This can be seen in Fig. |[ The nucleon binding 
is largest within the NL3 parametrization. In the TM1 case, the sl/2 + state is also deep, 
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Figure 3: Contour plot of nucleon densities for 8 Be without (left) and with (right) 
antiproton calculated with the parametrization NL3. The maximum density of normal 
8 Be is 0.20 fm -3 , while for nucleus with antiproton it is 0.61 fm -3 . 



but higher levels are less bound as compared to the NL3 and NL-Z2 calculations. This is 
a consequence of the smaller spatial extension of the potential in this case. The highest 
sl/2~ level is even less bound than for the system without an antiproton. 

For the antiproton levels, the TM1 parametrization predicts the deepest bound state 
with binding energy of about 1130 MeV. The NL3 calculation gives nearly the same 
binding, while in the NL-Z2 case, antiproton states are more shallow and have smaller 
spacing. It should be noted that the antiprotons are more strongly bound than was 
obtained in Ref. [|14]]. This follows from the fact that here we calculate both nucleons 
and the antinucleon self-consistently allowing the target nucleus to change its shape and 
structure due to the presence of the antiproton. The total binding energy of the system 
is predicted to be 828 MeV for NL-Z2, 1051 MeV for NL3, and 1159 MeV for TM1. For 
comparison, the binding energy of a normal 16 nucleus is 127.8, 128.7 and 130.3 MeV in 
the case of NL-Z2, NL3, and TM1, respectively. Due to this anomalous binding we call 
these systems Super Bound Nuclei (SBN). 

As a second example, we investigate the effect of a single antiproton inserted into the 
8 Be nucleus. In this calculation only the NL3 parametrization was used (the effect is 
similar for all three forces). The normal 8 Be nucleus is not spherical, exhibiting a clearly 
visible 2 a structure with the deformation (3 2 — 1-20 in the ground-state. Inserting the 
antiproton gives rise to compression and change of nuclear shape, resulting in a much less 
elongated nucleus with (3% ~ 0.23. Its maximum density increases by a factor of three 
from 0.20 fm -3 to 0.61 fm -3 . The cluster structure of ground state completely vanishes. 



A similar effect has been predicted in Ref. [18| for the case of the K bound state in 



the Be nucleus. In the considered case the binding energy increases from 52.9 MeV (the 
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experimental value is 56.5 MeV) to about 700 MeV! 



4 Life time, formation probability and signatures of 
SBNs 

The crucial question concerning a possible observation of the SBNs is their life time. 
The only decay channel for such states is the annihilation on surrounding nucleons. The 
mean life time of an antiproton in nucleonic matter of density ps can be estimated as 
t =< (JaVvc\Pb > _1 , where angular brackets denote averaging over the wave function of 
the antiproton and v re i is its relative velocity with respect to nucleons. In vacuum the 



NN annihilation cross section at low t> r ei can be parametrized as [21| a a = C + D/v Te \ 
with C=38 mb and D=35 mb. For < ps >— 2p this would lead to a very short life time, 
t ~ 0.7 fm/c (for v Te \ ~ 0.2). However, one should bear in mind that the annihilation 
process is very sensitive to the phase space available for decay products. For a bound 
nucleon and antinucleon the available energy is Q = 2m n — Bn — Bjf, where B N and Bjj 
are the corresponding binding energies. As follows from the calculations presented above 
(see Fig. this energy is strongly reduced compared to 2m N , namely, Q ~ 600 — 680 
MeV (TM1), 810-880 MeV (NL3) and 990-1050 MeV (NL-Z2) for the lowest antiproton 
states. 

For such low values of Q many important annihilation channels involving two heavy 
mesons (p, lu, rj, r)', ...) are simply closed. Other two-body channels such as irp, -kuj 
are considerably suppressed due to the closeness to the threshold. As is well known, the 
two-pion final states contribute only about 0.4% of the annihilation cross section. Even 



in vacuum all above mentioned channels contribute to a a not more than 15% | 22|| . There- 
fore, we expect that only multi-pion final states contribute significantly to antiproton 
annihilation in the SBN. But these channels are strongly suppressed due to the reduc- 
tion of the available phase space. Our calculations show that changing Q from 2 GeV 
to 1 GeV results in suppression factors 5, 40 and 1000 for the annihilation channels with 
3, 4 and 5 pions in the final state, respectively. Applying these suppression factors to 
the experimental branching ratios we come to the conclusion that in the SBNs the 
annihilation rates can be easily suppressed by factor of 20-30. There could be additional 
suppression factors of a structural origin which are difficult to estimate at present. This 
brings the SBN life time to the level of 15-20 fm/c which makes their experimental ob- 
servation feasible. The corresponding width, r ~ 10 MeV, is comparable to that of the 
cj-meson. 

Let us discuss now how these exotic nuclear states can be produced in the laboratory. 
We believe that the most direct way is to use antiproton beams of multi-GeV energy. This 
high energy is needed to suppress annihilation on the nuclear surface which dominates 
at low energies. To form a deeply bound state, the fast antiproton must transfer its 
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energy and momentum to one of the surrounding nucleons. This can be achieved through 
reactions of the type pp — > BB in the nucleus, 

p+(A,Z)^B + z(A-l,Z') , (1) 

where B = n,p, A, E. The fast antibaryon B can be used as a trigger of events where 
the antibaryon B is trapped in the nucleus. Of course, some additional soft pions can be 
emitted too. One can think even about producing an additional baryon-antibaryon pair 
and forming a nucleus with two antibaryons in the deeply bound states. In this case two 
fast nucleons will be knocked out from the nucleus. 

Without detailed transport calculations it is difficult to find the formation probability, 
W , of final nuclei with trapped antinucleons in these reactions. A rough estimate can be 
obtained by assuming that antiproton stopping is achieved in a single inelastic collision 
somewhere in the nuclear interior i.e. taking the penetration length of the order of the 
nuclear radius R . From the Poisson distribution in the number of collisions the probability 
of such an event is 

Wl = -exp^- y j, (2) 

where A^ 1 = a- in p Q and A -1 = (<7 in + cr^po (here <7 in and <ja are the inelastic and anni- 
hilation parts of the pN cross section). The exponential factor in Eq. (|2]) includes the 
probability to avoid annihilation. For initial antiproton momenta pi a b — 10 GeV we use 
a m ~ 25 mb, oa — 15 mb and get A ~ 1.6 fm which is comparable with the radii of 
light nuclei. For an oxygen target, using R ~ 3 fm leads to w\ ~ 0.17. 

In fact we need relatively small final antiproton momenta to overlap significantly with 
the momentum distribution of a bound state, namely, Ap ~ n/Rp, where Rp ~ 1.5 fm is 
characteristic size of antiproton spatial distribution (see Fig. []]). The probability of such a 
momentum loss can be estimated by the method of Refs. which was previously used 

for calculating proton spectra in high-energy pA collisions. At relativistic bombarding 
energies the differential cross sections of the pp — > pX and pp —> pX reactions are similar. 
The inelastic parts of these cross sections drop rapidly with transverse momentum, but 
they are practically flat as a function of longitudinal momentum of secondary particles. 
Thus, the probability of final antiproton momentum to fall in the interval Ap is simply 
Ap/pxab . For pi a b = 10 GeV and Ap = 0.4 GeV this gives 0.04. Assuming the geometrical 
fraction of central events ~ 20% we get the final estimate W ~ 0.17x0.04x0.2 ~ 1.4-10~ 3 . 
Even with extra factors ~ 0.1 which may come from the detailed calculations this is well 
within the modern experimental possibilities. 

Finally, we mention a few possible signatures of SBNs which can be used for their 
experimental detection. First of all, we remind the reader that according to the Dirac 
picture, any real antibaryon should be interpreted as a hole in the otherwise filled Dirac 
sea. Therefore, the nucleons from the positive-energy states of the Fermi sea can make 
direct transitions to the vacant negative-energy states of the Dirac sea. These super- 
transitions will be accompanied by the emission of a single pion or kaon depending on 
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the nature of the traped antibaryon. The energy of such a super-transition is fixed by the 
discrete levels of the initial and final baryons and according to our calculations should be 
of about 1 GeV. This 1-pion or 1-kaon annihilation is a unique feature of finite nuclear 
systems. In vacuum such transitions are forbidden by the energy-momentum conservation. 
Therefore, the observation of a line in the pion or kaon spectrum at energies between 1 and 
2 GeV would be a clear signal of the deep antibaryon states in nuclei. One can also look for 
narrow photon lines with energies in the range from 40 to 200 MeV corresponding to the 
transitions of nucleons and antibaryons between their respective levels. It is interesting 
to note that these signals will survive even if due to the lack of time the nucleus does not 
fully rearrange to a new structure. Another strong signal may come from the response 
of the nuclear remnant to the annihilation of the antibaryon in the deeply bound state. 
Since the remnant nucleus will initially be in a highly compressed state, it will expand 
and eventually break up into fragments. Therefore, the annihilation process will lead to 
rather cold multifragmentation with large collective flow of fragments. Both proposed 
signatures require rather ordinary measurements, which should be easy to perform with 
standard detectors. 

5 Discussion and conclusions 

Our main goal with this paper was to demonstrate that energetic antiproton beams can 
be used to study new interesting phenomena in nuclear physics. We discuss the possible 
existence of a completely new kind of strongly interacting systems where both the nucleons 
and the antinucleons coexist within the same volume and where annihilation is suppressed 
due to the reduction of the available phase space. Such systems are characterized by 
large binding energy and high nucleon density. Certainly, antinucleons can be replaced by 
antihyperons or even by antiquarks. We have presented the first self-consistent calculation 
of a finite nuclear system containing one antiproton in a deeply bound state. For this 
study we have used several versions of the RMF model which give excellent description of 
ordinary nuclei. The presence of an antiproton in a light nucleus like 8 Be or 16 changes 
drastically the whole structure of the nucleus leading to a much more dense and bound 
state. Even stronger effects are expected in the 4 He nucleus. It is clear however that these 
structural changes can occur only if the life time of the antibaryons in the nuclear interior 
is long enough. 

One should bear in mind that originally the RMF model was formulated within the 
Hartree and no-sea approximations. Implementing the Dirac sea may require serious 
revision of the model and inclusion of additional terms. Hartree calculations including 
the Dirac sea |14j] and Hartree-Fock calculations |l|, [19[] including exchange terms lead 
to smaller nucleon potentials in normal nuclei. Shallower potentials will produce smaller 
attraction for antinucleons, but the qualitative effect that the presence of antiprotons 
reduces repulsion and enhances attraction for nucleons will remain valid. We expect that 
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the additional binding and compression of the nucleus will appear even for an antinucleon 
potential as low as 200 MeV. 

Since nucleon densities in the considered systems could reach (2 — 3) p , it becomes 
questionable if the RMF model with nucleonic degrees of freedom is applicable at all. 
Coupling constants of RMF models, obtained by the fitting procedure, are predominantly 
constrained by observables at saturation density. On the other hand, the nuclear systems 
studied here are sensitive to the equation of state far above the saturation density. At 
such high densities nucleons might have to be substituted by quark degrees of freedom 
or modified by an admixture of them. Recently, we have used an extended version of the 
Nambu-Jona-Lasinio model to study bulk properties of systems composed of quarks and 
antiquarks |20[j . It has been found that deeply bound states also appear in this model, 



but the corresponding binding energies are considerably smaller by about a factor of 3. 

In summary, on the basis of the RMF model we have studied the structure of nuclear 
systems containing a few real antibaryons. We have demonstrated that the antibaryons 
act as strong attractors for the nucleons leading to enhanced binding and compression of 
the recipient nucleus. As our estimates show the life times of antibaryons in the nuclear 
environment could be significantly enhanced due to the reduction of the phase space 
available for annihilation. Narrow peaks in the pion or kaon spectra at the energy around 
1 GeV are proposed as the most clear signature of deeply-bound antibaryon states in 
nuclei. 
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